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Abstract We point out that the data on the onset of anomalous J/ip 
suppression as a function of nucleon numbers A and B could provide infor- 
mation on the dynamics of nuclear interactions. In particular the models of 
anomalous J /ip suppression by Blaizot and OUitrault (BO) and by Kharzeev, 
Lourenco, Nardi and Satz (KLNS) are based on different assumptions on the 
dynamics of nuclear collisions and lead to different predictions of the de- 
pendence of the onset of anomalous J/ip suppression on nucleon numbers 
of colliding nuclei. The data on this onset as function of A and B could 
become a tool for the study of the dynamics of nuclear collisions and bring 
further evidence on J/ ip suppression by new form of hadronic matter, possi- 
bly Quark-Gluon Plasma. In particular we propose to study J/ip suppression 
in A+Pb interaction with nucleon number of A between 58 and 73 or a bit 
higher. 

1 Introduction 

Experimental data of the NA38 and NA50 Collaborations at the CERN SPS 
[H, |, 1^, 1^ on J/'?/' suppression has shown that the anomalous J/ ip suppression 
is observed in near central Pb-|-Pb collisions but it is absent in peripheral 
Pb+Pb and in interactions induced by ions lighter (or equal) to Sulphur. 

Non- anomalous suppression of J/ip has been phenomenologically well 
described by the desintegration of J/ip by nucleons present in nuclei |[. 

There have appeared numerous attempts to describe the anomalous sup- 
pression of J/ip. We shall not try to make here a complete list of these 



1 



attempts. A review of approaches based on J/ip suppression by comovers 
can be found in Refs. [0, . For studies based on the initial state interactions 
see Refs.0 |ig. 

A rather abrupt increase of J/ip suppression with increasing transverse 
energy {Et) in Pb-Pb interactions observed by the NA50 Collaboration 
can also be interpreted as due to the formation of an "anomalous" state 
of hadronic matter, possibly of the QGP. Models of this type have been 
suggested by Blaizot and OUitrault (referred to below as BO) ||lT| and by 
Kharzeev, Lourengo, Nardi and Satz (KLNS) Both models have built 
in the onset of the formation of matter in the " anomalous" state so they are 
both able to describe the onset of the anomalous J / ip suppression at a given 
value of Et in Pb-Pb interactions at the CERN SPS. 

The purpose of this note is to point out that BO and KLNS models are 
based on different assumptions on dynamics of production of matter with 
high energy density in nuclear collisions and if they were the only possibili- 
ties it would be possible to learn which one gives the correct description by 
studying the occurence of anomalous J/ ip suppression as function of nucleon 
numbers of colliding nuclei. Since, in fact, BO and KLNS are not the only 
possibilities, the data on the dependence of anomalous J/ip suppression on 
nucleon numbers of colliding nuclei, if available, would provide an informa- 
tion on the mechanism of transverse energy production in nuclear collisions. 
In the Sect. 2 we shall discuss this point in more detail, in Sect. 3 we shall an- 
alyze the possibility of discriminating between the BO and KLNS models by 
using suitable experimental data and in Sect. 4 we shall present conclusions 
and comments. 



2 Nucleon number dependence of the onset 
of anomalous J/i/j production in BO and 
KLNS models 

In the present paper we shall use a simple model of nuclei as spheres with 
constant nucleon densities and radii of = 1, 2y4^/^fm where A is the nucleon 
number. In this simplified description some results can be obtained in a 
transparent analytical form. 

Both BO fill] and KLNS |jl2| assume that the mechanism of anomalous 



2 



J/ip suppression observed by the NA50 Collaboration at the CERN SPS is 
due to the dissolution of J /ip by the Quark- Gluon Plasma (QGP). The time 
spent by the J /if) in the QGP does not enter into these models what means 
that the dissolution is assumed to be immediate. The dynamics of at least 
the first decisive stage of the nuclear collision is assumed to be longitudinal 
in the sense that the criterion for the formation of QGP depends only on the 
particular tube- on- tube interaction. Following the picture of the Wounded 
Nucleon Model BO introduce density of participating nucleons per unit 



transverse area 

np{s,b) = TA{s)[l-exp{-aNTB{s-b))]+TBis-b)[l-exp{-crNTA{s))] (1) 

where ctat is the non-diffractive nucleon- nucleon cross- section, b is the impact 
parameter and s gives the transverse position of a nucleon with respect to 
the center of the nucleus A. Functions Ta{s) and Tb{s — b) give the nucleon 
densities per unit area in the transverse plane with 

/+00 
PA{s,z)dz (2) 
-oo 

The condition for the formation of QGP at a particular value of b and s is 
then stated in Ref.[|TT|l as 



[s, b) > nBo,c (3) 

where nBo,c is supposed to be just equal to the maximal possible value 
reached in S + U collisions. Note that the condition in Eq.(3) makes no 
explicit mention of the role possibly played by the formation times of the 
secondary hadrons or partons, which may be of some importance, see e.g. 
Ref.g. 



Neglecting exponential factors in Eq.(l), using the sharp sphere approxi- 
mation and multiplying np{s, b) by ctat =30mb to get a dimensionless quantity 
we obtain Eq.(3) in the following form 



i^BO = po'yN'ijRA - s2 + poaN2\/RB - {s- by > KbO,c (4) 



where kbo is defined as Up multipied by aN and po = 0.138/m ^. Following 



their earlier work |jT5[ BO assume in Ref.||TT| that all J/tp^s produced at 



values of b, s satisfying Eq. (4) are completely suppressed (independently of 
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whether they are produced directly or via e.g. radiative decays of x). When 
discussing the onset of QGP formation and thus also the onset of anomalous 
suppression of J/ip as function of nucleon numbers A and B we have to 
consider the collision of tubes with 6 = and s = 0, obtaining thus 

2poa7vl.2(Ai/3 ^ ^1/3) > ^^^^^ 

which in variables x = A^^^, y = B^^^ gives 

x + y = 1.006/tBo,c (6) 

In their model of the formation of QGP KLNS [0 use also the picture of 
tube- on- tube collisions, but in contradistinction to BO [|1T| they use as a 
criterion the interaction density [|16| defined (in the sharp sphere model) as 



i^K = ^ = ; / (7) 



The onset of J/ip suppression at the highest available Et (neglecting fluctu- 
ations in Et) corresponds again to 6 = 0, s = and Eq.(7) in this case leads 
to 

xy 



x + y 



what is equivalent to 



y = ha, a=1.00QKKc (8) 

X — a 

As pointed out by KLNS the kk as defined by the first part of Eq.(7) can 
in principle be determined experimentally. The number of nucleon- nucleon 
collision A^^ is proportional to the Drell- Yan (DY) pair production, whereas 
the number of secondary hadrons produced is proportional to the number of 
participating (wounded) nucleons. At the CERN SPS this statement is only 
approximately correct, since the data on transverse energy (Et) distributions 
and in particular the data on Et distributions associated to DY production 
do require some contribution to Et proportional to Nc, for details see the 
paper by Armesto, Capella, and Ferreiro [|| based on the Dual Parton Model 
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We would like to stress that the dynamics responsible for energy den- 
sity as given by Eq.(7) is rather different from that corresponding to Eq.(4). 
When considering the A+B interaction in the cms of individual nucleon- nu- 
cleon collisions, the numerator in Eq.(7) is proportional to the number of 
nucleon- nucleon collisions (within a given tube-on-tube collision) whereas 
the denominator is proportional to the sum of the lengths (and volumes) of 
the two interacting tubes. So one of possible interpretations of Eq.(7) is that 
Et relevant for the QGP formation is due to semihard gluon- gluon inter- 
actions, since the repetition of soft ones would be damped by the Landau- 
Pomeranchuk effect. The ratio in Eq.(7) is then proportional to the total 
transverse energy over the interaction volume and thus to energy density. 

For the sake of completeness we shall now present a few details concerning 
the calculations. When using the Gerschel and Hiifner 0, Q mechanism 
we have taken aabs =7mb and we have not introduced the effects due to 
the fluctuations of the Et = Exip) dependence. The calculations has been 
performed by using the expression 

S=Tr (9) 



Ra /-Stt 



iVn = / sds I de2xlRl - s^2Lf 
10 Jo 



where 



La{s) = - s2, Leih, s, 6) = ^Rl ~lP-s^ + 2bscos{9) 

when the expression under the square- root in is positive, and Lb(6, s, 0) 
when this expression is negative. 
The numerator is given as 

N= sds dee{Rl - - s^ + 2bscosie)) 
Jo Jo 

rLAis) fLBis,9) 

dzA / 

where we have put aa = 7mb. The last term is to be interpreted as B(k^^ — 
K^) or Q{k,c,bo — i^bo) respectively, or in the case of the two threshold scheme 
we have to make the replacement, see Ref. |jl2| 

Q{kc -k)^ OAQ{k\^ - kk) + 0.66(4, - ^^k) 



/ dzA / dzee-P^''^^^^^^^^'^ 

J-La(s) J-Ln{s,e) 
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where the first threshold corresponds to the dissolution of x by QGP and the 
second one to the dissolution of J/^. 

3 On the possibility of discriminating between 
BO and KLNS models by data 

In this section we shall show first that KLNS and BO schemes cannot be made 
equivalent in what concerns the dependence of anomalous J/V^ suppression 
on Et and on nucleon numbers A,B of colliding nuclei. As an illustration we 
shall consider a situation in which the onset of anomalous J /ip suppression is 
assumed to occur (without taking into account the fluctuations of Et) just 
in the interaction of the two longest tubes in S+U collision. In this academic 
example we take U-nucleus as spherical. In this case hbo,c = 9.312 and 
i^K,c = 2.086 and the curves for the onset of the anomalous J/ip suppression 
in the x = A^^^ and y = B^^^ plane are given by Eqs.(6,8) and this leads to 

(BO): y = 9.37-x, (KLNS): y = ^^ + 2.1 (10) 

The two curves are shown in Fig.l. As expected the two curves cross each 
other for values of x and y corresponding to S and U nuclei. 

With c=9.312 and kx,c=2.08 fixed by the assumption that anoma- 
lous J/ip suppression sets on just at the end of the E^ distribution for S-U 
interactions we can now study the onset of anomalous J/tp suppression in 
Pb-Pb collisions. In Fig. 2 we present J/ijj suppression calculated via Eq.(9) 
in Pb-Pb interactions. We have calculated ET{b) by using the relationship 
Erih) =0.325A^^(6)GeV, where N^ip) is the number of wounded (participat- 
ing) nucleons at a given value of h. 

Fig. 2 demonstrates an interesting, although very simple effect. When 
K^o c and kk,c are fixed by the condition related to the interaction of the two 
longest tubes in S-U interaction, the onset of the anomalous J/ ip suppression 
in Pb-Pb collisions appear at different values of Et in BO and KLNS schemes. 

Since, in fact, the anomalous J/^ suppression has not been observed in 
S-U interactions we find it more appropriate to determine the values kbo^c 
and nK,c by the data |Q on J/?/' suppression in Pb-Pb collisions. 

In Fig. 3 we present the NA50 Collaboration data @] on the survival prob- 
ability of J/ip a.s a function of the transverse energy Et in Pb-Pb interac- 
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tions at 158 GeV per nucleoli and three curves giving J/ ip suppression under 
different conditions: the curve denoted as "no anom." corresponds to pure 
Gerschel- Hiifner mechanism with (Jabs= 7mb, the curve "BO" is given by 
the BO scheme with kbo,c=9,75 and with all J/ip^s with kbo ^ i^bo,c com- 
pletely absorbed, the curve "KLNS" corresponds to the KLNS scheme with 
'^i<',c=2,43 and with all J/ip^s with kk > kk,c completely absorbed. 

In their work KLNS |]12| take into account that about 40% of J/ip^s 
observed in nuclear collisions is due to radiative x- decays and assume that x 
is suppressed at a different threshold in i^k than directly produced J/ip. This 
permits them to describe also the second decrease of survival probability of 
J/ip in the region of Et >100GeV indicated by the data. 

The same two threshold mechanism could also be implemented into the 
BO scheme and in this way both BO and KLNS schemes could describe the 
Pb+Pb data in a similar way. Within our simple model this could be achieved 
by increasing the value of aats and lowering the tresholds given by k,k,c and 
i^BO,c- We shall not follow this path here since the purpose of this note is to 
study only the onset (the first one) of the anomalous J/ip suppression and 
not the detailed shape of the survival probabability, including the possibility 
of a second threshold (the possible onset of x suppression). 

As can be seen in Fig.3 shapes of S{Et) obtained within BO and KLNS 
schemes in Pb-Pb interactions are rather similar when k,bo,c and kk,c are 
fixed by the condition that the onset of anomalous suppression starts at the 
same value of Et- It seems therefore that by using only the Pb-Pb data it 
would be very difficult to discriminate between the two schemes. 

In the rest of this section we shall try to find out whether one can dis- 
criminate between the BO and KLNS by studying the onset of anomalous 
J/ ip suppression as the function of atomic numbers A and B of the colliding 
nuclei. Note that we are neglecting the role of possible fluctuations in the 
energy density of matter formed in nuclear collisions, see Ref. ||18||. 

In Fig. 4 we show curves giving the onset of anomalous J/ip suppression 
in the x = A^^^ and y = B^^^ plane corresponding to values of c = 2.43 
and kbo,c = 9.75 (the values used also in calculations leading to Fig.3). 

As can be seen in Fig.4 for A-|-Pb collisions with 57 < A < 73 the survival 
probability as function of Et shows the presence of anomalous suppression 
within the BO scheme and no anomaly for the KLNS one. For nuclei with 
74 < A both KLNS and BO show the presence of anomaly but the anomaly 
sets on at different values of Et- As an illustration we show in Fig.5 the 



7 



results obtained again via Eq.(9) with kbo,c =9.75 and k,k,c =2.43 for the 
As^^-Pb interactions, still using Erib) =0.325GeViV^(6). 

4 Comments and conclusions 

As can be seen in Fig. 4 there do exist regions in the A,B plane where the 
BO picture of J/ip suppression based on the Wounded Nucleon Model []13| 



predicts anomalous J/ ip suppression whereas the KLNS scheme [|T2[ predicts 



its absence. For Pb as one of colliding nuclei this region for the other nucleus 
A contains nuclei with 

57 < A < 73 (11) 

including nuclei as Ni^Q, Cuq4 and Zuq^. For A > 74 the anomalous sup- 
pression is present in both schemes but sets on at lower value of Et in the 
BO case. 

To conclude: we propose to study experimentally anomalous J/i/j sup- 
pression at the CERN SPS in collisions of A-Pb or Pb-A with the nucleon 
number A within the interval given by Eq.(ll). The presence of anoma- 
lous J/ip suppression for A-Pb collisions with A within the interval given 
by Eq.(ll) would support the scheme of high energy density formation sug- 
gested by Blaizot and Ollitrault |Tl|, the absence of it would give evidence 



in favour of the Kharzeev, Lourengo, Nardi and Satz model ||T2 . 

A study of the onset of anomalous J/ip suppression in a series of A-Pb 
interactions including the available data on S-Pb and Pb-Pb, and some of 
A-Pb with A given by Eq.(ll) (or with A only a bit higher) at 158Gev per 
nucleon could bring interesting evidence on the dynamics of production of 
critical energy density for QGP formation in this energy range. Moreover, 
since at present the anomalous J/ip suppression has been observed in only 
one case, the second observation in A-Pb with a lighter A would support the 
idea of anomalous J/ ip suppression by an anomalous state of matter, possibly 
by QGP. As seen in Fig. 5 the difference between predictions of BO and KLNS 
schemes is not very large and some indeterminacy will be introduced by Et 



fluctuations [|12], |T8|, [T^ not considered here we still believe that a study of 
A-Pb interactions with A given by Eq.(ll) or nearby could bring in a very 
useful information. 
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Figure Captions 

Fig.l Curves denoted as KLNS and BO give lines of the onset of anoma- 
lous J/ip suppression in BO [|rT[] and KLNS |T2| schemes as functions of 



nucleon numbers A and B. The curves are given explicitly in Eq.(lO). The 
critical values are kbo,c=9.312 and c=2.086. 

Fig. 2 Dependence of the survival probability S{Et) on transverse energy 
Et calculated via Eq.(9) for KLNS and BO schemes with kbo,c=9.312 and 

c=2.086 (corresponding to the onset of anomalous suppression just above 
the S-U interactions) in the space of nucleon numbers A and B. The curve 
denoted as "no anom." gives the case with no anomalous suppression. 

Fig. 3 Comparison of survival probabilities of J/ip in Pb-Pb collisions in 
cases of no anomalous suppression ("no anom.") and BO and KLNS schemes. 
More details are given in the text. The data are taken from Ref. |Q and w.m.b 
means "with minimum bias". The critical values of k's are K5o,c=9.75 and 

Fig. 4 Curves giving the onset of anomalous J/ip suppression within 
KLNS and BO schemes with values kbo,c=9.75 and ^^-,£=2.43 obtained in 
comparison of KLNS and BO schemes with data in Fig.3. 

Fig. 5 Survival probability S{Et) calculated with kbo,c=Q-75 and /t/f^c=2.43 
for As''^-Pb interactions for BO and KLNS schemes. 
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